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I propose a quantum gravity model in which the fundamental degrees of freedom are information 
bits for both discrete space-time points and links connecting them. The Hamiltonian is a very 
simple network model consisting of a ferromagnetic Ising model for space-time vertices and an 
antiferromagnetic Ising model for the links. As a result of the frustration between these two terms, 
the ground state self-organizes as a new type of low-clustering graph with finite Hausdorff dimension 
4. The spectral dimension is lower than the Hausdorff dimension: it coincides with the Hausdorff 
dimension 4 at a first quantum phase transition corresponding to an IR. fixed point while at a second 
quantum phase transition describing small scales space-time dissolves into disordered information 
bits. The large-scale dimension 4 of the universe is related to the upper critical dimension 4 of 
the Ising model. At finite temperatures the universe graph emerges without big bang and without 
singularities from a ferromagnetic phase transition in which space-time itself forms out of a hot 
soup of information bits. When the temperature is lowered the universe graph unfolds and expands 
by lowering its connectivity, a mechanism I have called topological expansion. The model admits 
topological black hole excitations corresponding to graphs containing holes with no space-time inside 
and with “Schwarzschild-like” horizons with a lower spectral dimension. 

PACS numbers: 04.60.-m 


I. INTRODUCTION 

Three main approaches to the problem of quantum 
gravity have crystallized over the last thirty years. Super¬ 
string theory Jl| was introduced to bundle the graviton 
and gauge vector bosons in a single consistent mathe¬ 
matical framework. This, unfortunately requires many 
additional, as yet unobserved degrees of freedom as well 
as new symmetries and dimensions. It appears by now, 
however that strings themselves have to embedded in a 
larger framework, going under the name of M-theory . 
The major problem of M-theory, however is that it is not 
even properly defined at the non-perturbative level. 

Loop quantum gravity Q, instead, is a canonical ap¬ 
proach to the Einstein-Hilbert action formulated in terms 
of connections and vielbeins, analogously to the canoni¬ 
cal quantization of gauge theories. Correspondingly, the 
main variables are holonomies (loops). Despite its many 
successes and appealing structure at the Planck length 
fH, however, loop quantum gravity has not yet been able 
to successfully predict the geometries expected in the 
semi-classical limit. 

Finally, the Wilsonian approach based on causal dy¬ 
namical triangulations (CDT) [f| is the gravity equivalent 
of lattice gauge theories. It is formulated by discretizing 
space-time in terms of (generalized) triangles, formulat¬ 
ing the gravitational action via a modified Regge calcu¬ 
lus jb] and summing over all possible triangulations to 
obtain the quantum theory. The idea is to look for non- 
Gaussian ultraviolet fixed (UV) points that would define 
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quantum gravity in a non-perturbative sense, an idea go¬ 
ing under the name of asymptotic safety [7|. There is 
indeed compelling numerical evidence for the emergence 
of a de Sitter universe with finite Hausdorff dimension 
and scale-dependent spectral dimension [§]. The physics 
of these solutions is suggestive of a deep relation with a 
recent proposal |9j of a quantum gravity model in which 
time and space are fundamentally different at high en¬ 
ergies, with the Lorentz symmetry emerging only as a 
low-energy effective symmetry. This ties in with the fail¬ 
ure of an earlier approach based on Euclidean dynami¬ 
cal triangulations, with no imposed causal structure El- 
Euclidean dynamical triangulations have been recently 
revisited EH but no improvements upon CDTs have been 
found. 

Despite their differences, all these approaches share the 
common principle of starting from the correct classical 
Einstein-Hilbert gravity action and trying to reconcile it 
with quantum mechanics, either by enlarging the Hilbert 
space, by using different quantum variables or by dis¬ 
cretizing the degrees of freedom. In this paper I will 
advocate a different point of view by ascribing the non- 
renormalizability of perturbative quantum gravity to a 
complete misidentification of the fundamental degrees of 
freedom. I will suggest that geometric variables are not 
suited to formulate quantum gravity since the fundamen¬ 
tal degrees of freedom are purely combinatoric and de¬ 
scribe nothing else than information bits. In this picture, 
the geometry of our universe is an emergent property of 
the quantum ground state at a critical infrared (IR) fixed 
point while the small-scale physics is governed by an ul¬ 
traviolet fixed point for information bits with no geomet¬ 
ric space-time interpretation. Note that this is essentially 
the asymptotic safety scenario in disguise: in the usual 
approach one starts from the correct infrared fixed point 
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corresponding to general relativity and one looks for a 
non-trivial ultraviolet fixed point. Here I propose the 
reverse approach, to start from a different set of micro¬ 
scopic degrees of freedom describing the UV physics and 
trying to recover the existence of a flow to an infrared 
fixed point corresponding to general relativity. 

Of course I will not be able to complete this program 
in the present paper. However, I will try to substantiate 
my claim by introducing a surprisingly simple model of a 
network in which a 4-dimensional (4D) discrete universe 
is the ground state that emerges spontaneously in one 
of the three quantum phases of the model. This emerg¬ 
ing discrete universe is embodied by a low-clustering 12] 
graph with finite Hausdorff dimension 4. As in the CDT 
approach, the spectral dimension of this universe graph 
is smaller than its Hausdorff dimension. Admittedly, in 
this paper I concentrate only on the first step, from com¬ 
binatorics to topology, geometry and the link to general 
relativity and the Einstein-Hilbert action being still miss¬ 
ing. Nonetheless, I believe the phase structure of the 
emerging topology is very suggestive of a connection with 
quantized space-time. 


The model I propose is a generalization of Kazakov’s 
random lattice planar (2D) Ising model 0 to include 
dynamical topologies, governed by a second, antiferro¬ 
magnetic Ising model for edges of a random graph. In a 
nutshell, this second, edge Ising model renders the con¬ 
nectivity of the original Ising model dynamical by making 
it dependent on the dimensionless coupling constant. In 
Kazakov’s model the local limit at the critical Ising tem¬ 
perature describes 2D quantum gravity coupled to Majo- 
rana fermions, as can be derived exactly by matrix model 
techniques ]14{ ■ In the present model instead, the two rel¬ 
evant phase transitions do not depend on a critical tem¬ 
perature but are purely quantum. Here, the ground state 
of the original Ising model reacts to changes of the dy¬ 
namical connectivity 2d, which coincides with the dimen¬ 
sionless coupling constant of the model. Below the lower 
critical dimension di c = 1 the Ising model is disordered: 
correspondingly there is a first quantum phase transition 
at d = 1 in which the network breaks up into many small, 
disconnected and disordered components and any space- 
time interpretation is lost, the ’’universe” dissolving in 
fluctuating information bits. At the upper critical di¬ 
mension d uc = 4 the character of the Ising model ground 
state changes again. In this phase the Hausdorff dimen¬ 
sion of the ’’universe graph” coincides with its spectral 
dimension. Between these two quantum transitions there 
is a phase in which the spectral dimension is lower than 
the Hausdorff dimension, d s < du and the numerical ev¬ 
idence strongly points to a Hausdorff dimension taking 
the unique fixed value d h = 4, this being also supported 
by recent analytic results on the incommensurability be¬ 
tween the physical and the correlation length in finite size 
systems above the upper critical dimension fl5| . 

I cannot compute analytically the order of the two 
quantum phase transitions. However I conjecture that 
they correspond to UV (d= 1) and IR (d=4) critical 


points and realize thus the asymptotic safety scenario. 
The large scale dimension 4 of the universe would thus 
be predicted and related to the upper critical dimension 
of the Ising model. The UV critical point at d =1 would, 
additionally, imply that space-time ceases to exist at very 
small scales. As in the CDT approach [HJ the spectral di¬ 
mension of the universe decreases at small scales, but it 
does not stop at d=2, rather it goes all the way down to 
the Ising lower critical dimension d = 1 where space-time 
ceases to exist. 

The model can be treated at finite temperatures (T) 
in mean field theory. Of course, here one cannot think 
of temperature in terms of time averages of microstates 
since there is no time yet. Temperature, rather, has to be 
considered as a further external parameter/coupling con¬ 
stant that assigns non-vanishing probabilities to graphs 
different than the ground state according to their ener¬ 
gies, something like the stochastic temperature jig} as¬ 
sociated with the ParishWu fictitious time of stochastic 
quantization 0 - As I will show, this parameter can 
induce anisotropies in the universe graph, a most inter¬ 
esting possibility to distinguish a preferred direction. 

Mean field theory, however is unsuited to decide if the 
two quantum phase transitions are isolated or, rather, 
the endpoint of continuous lines of phase transitions at 
finite T. On the contrary, it should be quite accurate 
for high connectivities. Here it predicts a “topological 
cosmic expansion” due to a decrease in the graph con¬ 
nectivity when the temperature is lowered, i.e. the uni¬ 
verse graph actually “unfolds” in this phase rather than 
expand. There is no “big bang singularity”, the origin of 
the universe graph is a ferromagnetic phase transition in 
which the space-time information bits align and connect, 
i.e. discrete space-time itself self-organizes. Above the 
critical temperature there is only a hot information soup 
with maximal entropy and no possible interpretation as 
discrete space-time. During topological expansion the 
average distance on the graph scales exponentially with 
the temperature. A small drop in temperature causes a 
completely connected universe graph to unfold to a reg¬ 
ular graph of dimension 4. 

Finally, macroscopic “topological black holes” appear 
naturally as excitations of the model. It can be shown 
that such configurations are stable down to very low tem¬ 
peratures. 

The idea that in a proper theory of quantum grav¬ 
ity space-time should emerge spontaneously is not new. 
Models in which space is dynamically generated in quan¬ 
tum models with a well-defined notion of time were al¬ 
ready discussed in the 80s 0 . The same idea is also at 
the basis of the causal set approach to quantum gravity 
0 and of its recent development into energetic causal 
sets [go]. Ideas similar to those introduced in this paper 
have been previously considered in the so-called quantum 
graphity model [21]. The Hamiltonian of this model, 
however, is substantially different (and more complex) 
from the one proposed here and the derivation of the 
self-organization of a 4D universe-like graph from infor- 
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mation combinatorics is entirely missing. Another re¬ 
lated approach is based on structurally dynamic cellular 
networks (SDCNs) and is reviewed in in (22]. Finally, let 
me mention a recently developed, related approach [23} 
based on the idea of constructing discrete space-times by 
growing networks of simplicial complexes. 


II. THE MODEL 

Consider TV spin 1/2 information bits Si = ±1, for 
i = 1. . .TV and TV(TV — l)/2 spin 1/2 information bits 
Wij = 0,1, for i, j = 1... TV. A value of Si = +1 de¬ 
notes the existence of space time, while s* = — 1 indicates 
the absence of space time (or the presence of anti-space¬ 
time). A value = 1 denotes a connection between 
bits Si and Sj, a value Wij = 0 indicates that the two 
spins Si and Sj are not connected. 

All these information bits form a network with energy 

'Wik'Wkj 2 (i) 

k*. 

where the link variables are symmetric, Wij = Wji and 
vanish on the diagonal, w.a = 0 and where J is the unique 
dimensionless coupling. I use units in which K = 1, c = 1 
and all energies are measured in units of the standard 
Ising coupling, second term in ©, which is set to one for 
simplicity of presentation. 

The second term in this energy function is the stan¬ 
dard ferromagnetic Ising model. If the links Wij would 
be uniformly drawn from random adjacency matrices of 
degree 4, the model would be Kazakov’s random lattice 
Ising model in two dimensions [13] . The first term in the 
energy functional, on the other side, is simply a nearest- 
neighbours (sharing a common vertex) antiferromagnetic 
Ising model for the link spins. The generalizations with 
respect to Kazakov’s model, thus consist in dropping the 
restriction to degree 4 and drawing the random adjacency 
matrices from a Gaussian distribution. 

The competition between the vertex ferromagnetic 
coupling and the link antiferromagentic one creates ” link 
frustration” in the model. Indeed the vertex ferromag¬ 
netic coupling favours the creation of many links (posi¬ 
tive values of aiq) in a state with the majority of vertex 
spins aligned, corresponding to an incipient space-time: 
space-time points have a tendency to link together. How¬ 
ever, due to the antiferromagnetic link coupling, creating 
many links costs energy and is energetically disfavoured. 
It is particularly nearest neighbour links, sharing a com¬ 
mon vertex, that are disfavoured. As I now show, the 
compromise is to create links between space-time points 
but to avoid triangles involving nearest neighbours links. 
The result is a graph with power-law extension, exactly 
what one would expect for a discretized universe. The 
dimension of this universe, encoded in the power-law for 
the average distance is determined by the unique dimen¬ 
sionless coupling J of the model. 


To show this, let me consider a /c-regular graph, i.e. a 
configuration with Sj = +1, Vz and such that each vertex 
has exactly k incident edges (degree k in graph parlance 
0). When k is even, such a graph resembles locally 
a fc/2-dimensional lattice. I will now show that such a 
graph is a local minimum of the energy ©• 

To this end, suppose we are given such a graph and let 
me first change one single vertex spin from s» = +1 to 
Si = — 1. The corresponding energy change AEi = k > 0 
is positive. Changing vertex spins, thus costs energy. 
Let us now try to eliminate an existing connection, i.e. 
changing = 1 to Wij =0. In this case there are con¬ 
tributions from both terms in the energy function. Be¬ 
fore the elimination, the existing connection contributed 
Eij = J(k — 1) — (1/2). After the elimination, of course 
the contribution of this connection vanishes. The en¬ 
ergy change due to the elimination of a connection is 
thus A e ii m Eij = (1/2) — J[k — 1). Let me now consider 
adding a previously non-existent connection Wij . In this 
case it is the energy contribution before the addition that 
vanishes, while, after the addition I have an additional 
energy E^ = J/c — (1/2). The energy change for adding a 
connection is thus A ac jd Eij = Jk — (1/2). Requiring that 
both eliminating and adding a connection costs energy 
we obtain the stability condition 


1 

2 J 


< k < 1 + 


2 J 


( 2 ) 


To proceed, let me compute the total energy of a k- 
regular graph on TV vertices. This is easily obtained as 


ATv.fc = N 


J 


k[k - 1) 



(3) 


This expression is minimized when the vertex degree 
takes the value k = 1/2+ (1/2J) and this value of k 
satisfies the stability condition ©. Defining 


J ~ 4d-l 


(4) 


I have obtained the result that for any choice of integer d 
a 2d-regular graph is a local minimum of the energy and 
that this minimum is the one of lowest energy among 
all possible regular graphs. As anticipated, the links ar¬ 
range themselves to form a locally lattice-like graph with 
2d edges at each vertex; the number of triangles is min¬ 
imized. Thus, d plays the role of the ’’bare” dimension 
of the universe. I call this dimension ’’bare” since, as I 
will show below, the physical Hausdorff dimension of the 
universe can differ from d. 

Of course, I have not shown that 2d-regular graphs 
are the true global minima of the energy functional ©■ 
There could be other, non-regular graphs, with even 
lower energies. That this is not so can be made plau¬ 
sible with the numerical methods introduced in the next 
section. 
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III. THE PROPERTIES OF EMERGENT 
UNIVERSE GRAPHS 


The energy function © defines a Maxwell-Boltzmann 
distribution on the space of graphs. This distribution 
is peaked around the ground state of H , which corre¬ 
sponds to the semi-classical “universe graph”. Graphs 
with large weights in the distribution correspond to fluc¬ 
tuations about this semi-classical universe graph. In the 
rest of this paper I shall focus on the properties of the 
semi-classical universe graph as a function of the dimen¬ 
sionless coupling J. 

To do so I will reverse a technique used to analyze 
neural networks such as the Hopfield model [25j]. There, 
a network dynamics is posited and an energy functional 
is looked for that acts as a Lyapunov function, i.e. such 
that the energy minimum determines the fixed points 
of the network dynamics. Here I will do the contrary, 
namely I will look for a network dynamics whose fixed 
points correspond to minima of the energy ©■ 

To this end let me start form a random initial configu¬ 
ration and sequentially update the vertex and link spins 
according to the rule 


Si(t + 1) = sign + (hi{t)) , 

h i = J2w ij s j , (5) 

Wij(t + 1) = 0± ( hij ) 

hij = i SiSj - W - T Wik ^ ’ ( 6 ) 

k^j k^i 

k^i k^j 


where 0 denotes the Heaviside function and ha = 0. The 
subscripts ” +” and ” ±” on the sign and 0 functions indi¬ 
cate what is the rule to follow in the (rare) cases in which 
the argument is zero. Complete symmetry would require 
a totally random choice in these cases. This, however 
can lead to a two-component network in which space- 
time vertices have the opposite sign in the two compo¬ 
nents. Since I am interested in the emergent properties of 
just one of these components I will maximize its size by 
making an asymmetric choice for the sign function deter¬ 
mining the orientation of the space-time vertices. This 
is embodied in the subscript ”+”, which indicates that, 
in this case, the + sign has to be chosen. This favours 
positive vertex spin values so that one-component graphs 
are preferred. The subscript ”±” on the Heaviside func¬ 
tion, instead indicates that a purely random choice will 
be made when the argument vanishes. This choice also 
helps to avoid being stuck in sub-optimal minima during 
the evolution. 


Ei to the energy changes according to 

Ei(t + 1 ) = —Si (t + 1 ) ^2 u ’ij ( t ) s j (*) 

= —sign + ( hi(t )) hi(t) = -\hi(t)\ 

< —Si(t)hi(t) = Ei(t) . (7) 

With the exact same procedure for the update of a link 
spin hij we obtain 

Eij(t+ 1 ) = -Wij(t+l)hij{t) = - 0 ± ( hij(t )) hij(t) . ( 8 ) 


Now we have two possibilities 

• hij(t) < 0 -a Eij(t + 1 ) = 0 < - Wij(t)hij(t ) = 
Eij (t) 

• hij(t) > 0 -a Eij(t + 1) = - hij(t ) = ~\hij(t)\ < 
-Wij(t)hij(t) = Ejj (t) 

and in both cases we obtain Eij(t + 1) < Eij(t ), which 
proves the claim that the energy cannot increase during 
sequential evolution (note that, for both the sign and the 
Heaviside functions, the exact procedure on how to treat 
the undefined cases in which the argument vanishes has 
no effect on this result). As a consequence, every mini¬ 
mum of the energy © is a fixed point of the sequential 
network evolution ([5]) and ©. 

Every such fixed point defines an (undirected) graph. 
Several classes of graphs have been studied [E|] exten¬ 
sively. The simplest such one is the class of regular 
graphs, for which every node has a fixed degree fc, i.e. 
a fixed number k of neighbours. Regular graphs are 
typically highly clustered, with the notable exception of 
d-dimensional hypercubic lattices, which have vanishing 
clustering coefficient since they do not contain any trian¬ 
gles, and in which the average distance scales as N 1 ^. 
This is the exact opposite of classical random graphs [26} 
in which connections between nodes are chosen randomly 
with a given fixed probability. Random graphs have small 
clustering c rg = ( k)/N where (k) is the degree expecta¬ 
tion value and small average distance that scales loga¬ 
rithmically with N (the so-called small world effect). In 
classical random graphs the degrees follow a Poisson dis¬ 
tribution. In order to obtain graphs that better model 
real-world networks a class of generalized random graphs 
was introduced [27}, in which connections are chosen ran¬ 
domly subject to the constraint that the degree distribu¬ 
tion must follow a predetermined law. In generalized 
random graphs the clustering depends, of course, from 
the chosen input degree distribution. Actually, only the 
first two moments of this distribution matter, 


1 «fc 2 )-(fc) 2 ) 2 

Cgrg “ N (k } 3 


(9) 


I will now show that the energy function © cannot 
increase along a sequential evolution (0 and ®: E(t + 
1) < E(t). To this end let me begin by considering the 
update of vertex spin s,;. The corresponding contribution 


This reduces to c rg for a Poisson distribution, as ex¬ 
pected. Generalized random graphs are also small 
worlds, the average distance remains a logarithmic func¬ 
tion of N. These two classes describe static graphs, with 
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a fixed number N of nodes. The number of nodes grows, 
instead in the class of evolving grap hs f28[ . Many evolu¬ 
tion rules have been considered ©}, essentially all aimed 
at generating scale-free graphs with power-law degree dis¬ 
tributions of various forms. 

The graphs that minimize the energy © constitute yet 
another, novel category for which 


which can be also written as c = number of closed triples 
/ number of connected triples. In this last form it can 
be easily evaluated once the adjaceny matrix A of the 
ground state graph is known: 


Tr A 3 

l' A 2 ' 


( 10 ) 


• The number of nodes is not growing but fixed. 

• There is no random assignment of connections but, 
rather, these are dynamically determined. The 
dynamics, however involves both connections and 
nodes. 


The spectral dimension d s measures the connectivity 
of the graph |29[, the dimension that a particle moving 
on the graph would feel. It is defined via the scaling of 
the return probability p r {t) to the initial point after t 
steps of a random walk on the graph [30]. 


• The dynamics is not an ad-hoc rule of how a graph 
is grown but is governed by an energy minimization 
principle. 

• Contrary to generalized random graphs and typical 
evolving graphs, the degree distribution is not an 
input but an output. 

I shall call such graphs dynamical graphs and, in the 
following, I shall describe their main properties. 

In order to investigate these dynamical graphs, I have 
repeatedly started with random values of all the spins 
and let them evolve according to the dynamics © and © 
until a fixed point is reached. The properties of interest 
can then be extracted from the adjacency matrix of this 
graph. Finally, an average over the results form different 
runs is taken. Of course, care has to be taken to avoid 
ending up in fixed points corresponding to metastable, 
local minima of the energy function. Fortunately, these 
typically involve configurations for which not all space- 
time information bits take the same value +1 and can 
thus be recognized easily. Actually, it could also be possi¬ 
ble that a fixed point cannot be reached since the ground 
state is degenerate and so the network dynamics would 
wander among the different ground states with the same 
energy. This would be no problem however, since any¬ 
one of these ground states would correspond to a viable 
universe graph. 

As expected, I have observed that, for d integer, the 
minimum energy configuration always corresponds to a 
2d-regular graph with exactly 2d edges for every vertex. 
By the degree sum formula 2e = JA> 3 i Vi, with e the 
number of edges and Vi the number of vertices of de¬ 
gree i, we can derive that these graphs have exactly dN 
edges. The connectivity details of the ground state graph, 
like clustering coefficient, number of 4-cycles and others 
can indeed vary from run to run but they do correspond 
to finite size effects: the bulk properties convergence to 
unique values for large values of N. 

The three most important quantities that I have mea¬ 
sured are the clustering coefficient , the spectral dimension 
and the Hausdorff dimension of the emergent dynamical 
graphs. As already anticipated, the clustering coefficient 
c ©I is related to the frequency of triangles in the graph: 
c = 3 x number of triangles / number of connected triples, 


p r (t) ~ t da/2 . (11) 

For infinite graphs this scaling relation is valid in the 
limit t —» oo. For finite ^-regular graphs, the return 
time to the initial point is t = N ©j and thus p r {t) —> 
1 in the limit t oo. The correct scaling is typically 
found in the intermediate region 1 <f < 0(N 4: ^ k ) where 
finite size effects are suppressed [3(ij |. To compute the 
spectral dimension one introduces the degree distribution 
di denoting the number of edges at each vertex i. The 
matrix 


Mij — d A.ij , (12) 

represents the transition probabilities from vertex i to 
vertex j in a random walk on the ground state graph. 
The average over all possible initial points of the return 
probability of the random walk after t steps is then given 

by 


Pr(t) = ^ TrM * ■ 


(13) 


The best fit of this to a power t~ ds / 2 over the appropriate 
f-interval defines then the spectral dimension d s . 

Finally, the intrinsic Hausdorff dimension dn measures 
how t he g raph volume scales with distances T> on the 
graph ©I, 


V ~ N 1,dH . 


(14) 


Of course, this depends on the exact definition of the 
graph distance. Here I will follow 3l| and take T> = (D) 
as the average distance among vertices of the graph. This 
can be also computed from the adjacency matrix A of the 
ground state graph. Indeed, the minimal graph distance 
between two vertices i and j is given by the smallest 
integer power t-ij such that (A i ’'>) .. ^ 0. The average 
distance on the graph is thus given by 


(D) = 


N(N - 1) 




(15) 


Let me stress that this is an intrinsic property of the 
graph that does not need an embedding in an extrinsic 
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FIG. 1: The average number of disconnected components of 
dynamical space-time graphs as a function of the coupling 
constant d for N = 120. All graphs with values between 
d = 0.75 and d = 1.25 have connectivity 2, corresponding 
to a one-dimensional Ising model for the space-time vertices. 
This supports the existence of a quantum phase transition at 
the lower critical dimension d = 1 of the Ising model. 


Euclidean space to be defined and measured. It is the 
dimension that would be measured with clocks and rods 
by ’’inhabitants” of the graph. This dimension is simply 
extracted as the best fit of (fl5l) to a power of TV as in 

m- 

For large values of N the convergence to the ground 
states slows down considerably and so does the comput¬ 
ing time needed to evaluate the quantities of interest. 
This difficulty is compounded by the increasing probabil¬ 
ity of ending up in local minima of the energy function. 
These difficulties, combined with the limits of the avail¬ 
able computational resources have bounded the accessi¬ 
ble size of graphs to vertex numbers of up to a maximum 
of TV = 250 (TV(TV - l)/2 = 31125 possible links). 

The first important observation is the strong evidence 
for a first quantum phase transition at d = 1 (J = 1/3), 
a quantitative change in behaviour of the ground state 
at zero temperature as a function of the model coupling 
constant [321 ]. When d = 1 is approached from above the 
dynamical graphs start to break up into many discon¬ 
nected components, as shown in Fig. 1 for TV = 120. 

This happens from d = 1.25, where some vertices of 
the graph start to have connectivity 2 (corresponding 
to one dimension). Of course, the number of connected 
components diverges like TV for d < 0.75, where some 
vertices begin to have connectivity 1 (i.e. less than one 
dimension). The interesting region, however, lies be¬ 
tween d = 1.25 and d = 0.75, where the connectivity 
is always 2 (corresponding to the one-dimensional Ising 
model) for all vertices and the TV-dependence is much less 
pronounced. Here, there is a clear evidence of a quantum 
phase transition at d = 1, where the number of discon¬ 
nected components increases dramatically even at these 
small finite sizes. Note that this happens even if the ficti¬ 
tious energy-descent dynamics is chosen to favour aligned 
space-time spins, as in eq. ©■ For a completely symmet¬ 


ric dynamics the proliferating disconnected components 
have also disordered space-time spins. 

In order to understand the nature of this phase transi¬ 
tion it is helpful to consider the model © as a standard 
nearest-neighbours Ising model with dynamical dimen¬ 
sionality d. As is well known [24}, the Ising model has 
a lower critical dimension d\ c = 1 below which it is al¬ 
ways disordered (at zero temperature). At this lower 
critical dimension domain walls lower the energy and it 
becomes thus energetically favourable for the system to 
create many regions of spins of different signs, eventually 
completely disordering the system. In the present model 
with dynamical edges, however, there is an alternative 
way to lower the energy, namely simply breaking up the 
links instead of creating domain walls. This is exactly 
what is happening. This first quantum phase transition 
at d = 1 is thus associated with the lower critical dimen¬ 
sion of the Ising model. 

Dynamical graphs, emerging as ground states of ©, 
are low-clustering graphs. Their clustering coefficients 
scale as 1/TV, as is the case for random graphs and gen¬ 
eralized random graphs. Any additional power posited in 
the TV-dependence resulted in statistically non-significant 
parameters of the non-linear regression. In the following 
table, the clustering coefficients Cd g of dynamical graphs 
and their regression (as a function of N) standard devi¬ 
ations are reported together with the clustering coeffi¬ 
cients of generalized random graphs with the same, uni¬ 
form degree distribution, obtained from ©. 


d 

2 

3 

4 

5 

NCdg 

2.4 

4.52 

6.27 

8.39 

& Nc 

0.15 

0.09 

0.22 

0.17 

NCg r g 

2.25 

4.17 

6.125 

8.1 


The two values are very close, although the values for 
dynamical graphs are slightly higher. This may be the 
result of the limited number of vertices in the numerical 
simulations. In any case it is evident that clustering is a 
finite size effect in dynamical graphs. 

In Fig. 2 1 plot together the spectral and Hausdorff 
dimensions as a function of the bare dimension d, the 
coupling constant of the model. As expected, the spectral 
dimension d s coincides with the bare dimension d. Not 
so for the Hausdorff dimension du though. 

Indeed, Fig. 2 highlights the evidence for a second 
quantum phase transition at d = 4 (J = 1/15). Al¬ 
though the results are slightly higher than expected, the 
indication is that, for d > 4 Hausdorff, spectral and bare 
dimension coincide while, in the range 1 < d < 4 the 
spectral dimension is lower than the Hausdorff dimen¬ 
sion, which is fixed at the value dn = 4, the required 
value to describe our observed space-time. The differ¬ 
ence dn — d s can then be considered as a geometric dis¬ 
order parameter for this quantum phase transition and 
is shown in Fig. 3. 

There is actually a plausible explanation of the physi¬ 
cal origin of the special role of the dimension 4. To this 
end let me consider again the interpretation of model © 
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FIG. 2: spectral and Hausdorff dimensions as a function of 
the coupling constant d, showing the decoupling of these two 
dimensions at d = 4. 



FIG. 3: the geometric disorder parameter dH-d s as a function 
of the coupling constant d. 


as a standard nearest-neighbours Ising model with dy¬ 
namical dimensionality d. The crucial point is that the 
Ising model has also an upper critical dimension d uc be¬ 
yond which fluctuations can be neglected altogether and 
mean field theory becomes exact [24]. For the Ising model 
this upper critical dimension turns out to be d uc = 4. 

It is often stated that hyperscaling fails above the up¬ 
per critical dimension. This is because, in the standard 
approach to finite size scaling [33, HH , one considers the 
correlation length as bounded by the system size L. It 
has been shown fl5j . however that hyperscaling can be 
fully restored if this assumption that the finite size corre¬ 
lation length be bounded by the length scale is relaxed. 
The correct scaling is 


£ oc L d/d " c , d > d uc , 

£ oc L , d < d uc . (16) 


This is because, in finite size samples, the correct correla¬ 
tion length is not commensurate with the physical length 
above the upper critical dimension [TrJ . It coincides with 
the physical length below the upper critical dimension 
but differs above it. One can now define also two volumes, 
V = L d and the effective volume V e g = £ d “°. These, in¬ 


stead, coincide above the upper critical dimension but 
differ below it. An alternative, but equivalent way to 
represent the incommensurability found in [l5[ when the 
number of spins, rather than the geometric size, is fixed, 
is to consider the effective volume as the volume actually 
dynamically occupied by the system in any number of di¬ 
mensions. In terms of this effective volume, the physical 
length scales as: 

L oc Veff d , d > 4 , 

L oc V^ff 4 , d<4, (17) 

where I have already explicitly used d uc = 4 for the Ising 
model. This is exactly what I have found for the Haus¬ 
dorff dimension of model ©• In other words, the quan¬ 
tum phase transition of © at d = 4 (or equivalently at 
J = 1/15) is due to the change in critical behaviour of the 
Ising model at its upper critical dimension. For d > 4 the 
model is thus in a mean field phase with d = d s = dn, 
for d < 4, instead fluctuations become important and 
the model is in a fluctuations dominated phase in which 
dn = 4 and d s < dn- The dimension 4 of the universe 
would thus be simply related to the fact that 4 is the 
upper critical dimension of the Ising model. Note that a 
spectral dimension lower than the Hausdorff dimension 
is also exactly what is predicted by the CDT approach 
[§:]. On the graph sizes accessible to numerical simula¬ 
tions, I have not been able to observe a scale-dependence 
of the spectral dimension, as predicted by the CDT ap¬ 
proach [8|. It seems, however natural to conjecture that 
the quantum phase transition at d = 4 corresponds to an 
infrared-stable critical point, especially in view of the fact 
that the Gaussian model is an infrared fixed point of the 
Ising model at d = 4. In this case the model © would 
predict the large-scale dimension 4 of the universe and re¬ 
alize at the same time the asymptotic safety scenario of 
quantum gravity Q. The other quantum phase transition 
at d = 1, in fact, would necessarily imply the existence 
of an ultraviolet-stable fixed point either at d = 1 or in 
the intermediate region 1 < d < 4. In both cases the 
spectral dimension would be scale-dependent. The latter 
case would correspond to a universe of lower dimensional¬ 
ity at short distances. The former, even more interesting 
case, would imply that there is no space-time at all at 
short distances, only fluctuating information bits, some¬ 
thing like the “quarks of quantum gravity”. The fact that 
T = 0 is an ultraviolet fixed point of the one-dimensional 
Ising model f36j supports this picture. 

To conclude this section I would like to discuss the 
role of quantum mechanics in the present model. Indeed, 
the Hamiltonian © appears at first sight as a purely 
classical one. But so is also the Hamiltonian of the orig¬ 
inal 2 D Kazakov random lattice Ising model jl3j. Kaza¬ 
kov’s seminal insight was that a new quantum degree of 
freedom, the fluctuation of the metric, i.e. 2D quan¬ 
tum gravity, emerges at the critical point in the local 
limit. Indeed, it is known since the old days of stochastic 
quantization that the Euclidean quantum par¬ 

tition function of a field theory can be obtained as the 








thermal partition function of the corresponding classical 
system coupled to a fictitious thermal bath. The main 
idea of the present model is essentially the same. The 
difference lies in the fact that, here, the critical points 
for which I have given first evidence are two and both 
lie at zero temperature. When taking the T —> 0 limit 
of the partition function corresponding to ©, the con¬ 
figuration sum reduces to a sum over the many vertex 
and graph configurations of the same minimal energy. 
If no transverse magnetic fields are applied, model (HD 
is identical at the classical and quantum levels since all 
cr 3 Pauli matrices attached to vertices and edges com¬ 
mute. The graph configurations of minimal energy that 
survive in the T = 0 partition function can thus be taken 
to represent the residual zero-temperature quantum fluc¬ 
tuations. The idea is that 4 D quantum gravity emerges 
when these become critical in the IR at d = 4. The other, 
UV critical point at d = 1 would then represent the quan¬ 
tum behaviour of the theory at very small scales, where 
space-time dissolves into random information bits. Note 
that the Hamiltonian © represents one of the typical 
paradigms for quantum phase transitions [32|. Typically, 
a Hamiltonian H = Hq + JHi can display non-analiticity 
in its ground state if Hq and Hi commute, exactly as in 
the present model, leading to level-crossings at particu¬ 
lar values of J while the wave functions remain the same 



IV. TOPOLOGICAL EXPANSION 

In this section I will focus on the finite-temperature 
behaviour of the model ©■ As stressed in the intro¬ 
duction, here temperature cannot be thought of in terms 
of time averages of microstates, since there is no time 
yet. Rather, it has to be considered as a further exter¬ 
nal parameter that assigns non-vanishing probabilities to 
graphs different than the ground state according to their 
energies. 

I will adopt again a technique that is used in the anal¬ 
ysis of neural networks [25| . This consists in promot¬ 
ing the deterministic network dynamics © and © to a 
stochastic rule for the probabilities of assuming the bi¬ 
nary values, 

Prob (si(t + 1) = +1) = / ( hi(t )) , 

Prob + 1) = +1) = / (hij(t)) , (18) 

where 

fW = ! + 1-20H > ( 19 ) 

is the Fermi function at temperature T = 1//3. I will 
now show that the equilibrium reached by this stochastic 
update rule is the Boltzman distribution corresponding 
to the energy function ©• Therefore, assuming an equi¬ 
librium configuration, the two above rules will provide 
coupled mean field equations for mean space-time ’’mag¬ 
netization” < s > and the average degree < k > of the 


graph vertices. This defines a temperature-dependent ef¬ 
fective dimension of the universe graph. I don’t expect 
mean field theory to be accurate enough to derive if the 
two quantum phase transitions are isolated or, rather, 
the endpoints of continuous phase transitions as a func¬ 
tion of the temperature. However, it should be accurate 
enough for high connectivities < k >. 

Let me denote by F s ( Si -A — s,) the probability of a 
vertex spin flip and by F w (Wij —> 1 — lOy) the probabil¬ 
ity of a link spin flip. As a consequence of (1T5T) these are 
given by 


F s ( Si -A —Si) 


Fyj {wij ? 1 rCij) 


e ~/3hiSi 
2 cosh (/3/iiSi) 

g—fihijSij 

2 cosh (phijSij) ’ ^ 


where I have introduced Sij = 2 (w^ — 1/2). Let me de¬ 
note by D (... Si... ...) the probability distribution 

for activity patterns of all vertex and link spins. Due to 
the vanishing of the diagonal terms, wa = 0, and sym¬ 
metry = Wji of the link variables I will consider D as 
function of w^ for i < j only. In equilibrium (denoted 
by ”e”), the distribution D e must satisfy the conditions 


F s (Si -A -Si) D e (.. .Si... Wij ...) 

= F s (-Si -A Si) D e (- Si...Wij ...) , 

Fw (Wij = 0 -a Wij = 1) D e (... Si ... Wij = 0... ) 

— Fyj (Wij — 1 ^ — 0) D e (. . . Si . . . Wij — 1 . . /^l.) 

Using (®D one obtains readily 

D e (■ • ■ Sj . . ■ Wij . . ■ ) _ ^2phiSi 

D e (• • * Si .. . Wij . . . ) 

D e (... Sj... w^ = 1...) _ phn ^22) 

D e ( - Sj.. .w^ = 0 ...) 

Finally, one easily verifies that these conditions are met 
by the equilibrium probability distribution 

l)r(...S,...W,j...)- 1 f S// 5 (23 ) 

where H is the energy function © and Z the correspond¬ 
ing partition function. The factor of 2 in the exponent of 
(l22l) arises since all fundamental variables, both vertex 
spins Si and link spins Wi<j appear twice in the energy 
function, the latter due to the symmetry w^ = Wji. The 
equilibrium distribution is thus the Boltzmann distribu¬ 
tion corresponding to the energy function ©. 

Having established that m are the probabilities in 
the equilibrium Boltzmann distribution at inverse tem¬ 
perature /? I can use them to compute the expectation 
values of both vertex and link spins in thermal equilib¬ 
rium. Let me begin by the vertex spins. In this case, the 
expectation value is 


(Si) 


+1/ (hi)-lf (-hi) = tanh 



• (24) 
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In the mean field approximation, which becomes exact 
in large number of dimensions, I can bring averages (...) 
inside non-linear functions. With the usual homogeneity 
Ansatz ( Si) = (s), independent of i, I obtain 

(s) = tanh (/3(fc)(s)) , (25) 

where I have introduced the mean degree (ki) at vertex 
b 


(fo)=^(wy), (26) 


and assumed that it is independent of i, (ki) = (k), Vi. 

One proceeds similarly for link spins. In the mean field 
approximation, the thermal average is given by 

(wij) = f(hij) = f Q(s) 2 - J{k) + J(wij)j • (27) 


In this case, however, I am not interested in homogeneous 
solutions in the link variables, but, rather, only the av¬ 
erage degree (E7>1) is expected to be independent of the 
vertex label. I will therefore sum the left hand side of 
(EZl) over all j ^ i to obtain 



At high temperatures one can expect the thermal average 
of the link spins to approach its limiting value 1/2 from 
below. I will thus approximate inside the Fermi function 
(wij) ~ 1/2 — eij, Cij > 0 and neglect the deviations e^, 
keeping however in mind that they are positive. This 
gives 


(fc) 

N- 1 


r 0 (s) 2 -j(V) 



(29) 


where the superscript means that 0 must be ap¬ 
proached from below in the Fermi function. Equations 
(l25l) and (l29l) constitute the coupled mean field equations 
for ’’space-time magnetization” and mean degree, respec¬ 
tively. They describe topological fluctuations of discrete 
space-time graphs and their connectivity at inverse tem¬ 
perature j3 and should provide a correct picture of these 
emerging universe graphs as long as (k) 1. 

To solve these equations I will introduce the new vari¬ 
able x = (1/2)(s) 2 — J ((k) — (1/2)) and the new linear 
function 

!'W = ]v^T(^W 2 + 5)-7(FrT) 1 - < 30 > 

and rewrite the mean field equations as 


(s) =tanh (f3(k)(s)) , (31) 

g(x) =f~(x). (32) 

When f3(k) < 1, eq. (15TT) has the unique solution (s) = 0. 
This is the case in the high temperature regime /3 —>■ 0 


(for finite ( k )). In this limit, the Fermi function in 
(15211 reduces to its limiting value 1/2, giving the solution 
g(x) = 1/2 or (k) = (k ) ma x = (N — l)/2, which is indeed 
finite. This is the information phase of the model, a hot 
soup of information bits with no space-time interpreta¬ 
tion. When P(k) > 1, instead, two solutions of (1511) exist, 
the same solution (s) = 0 as before and a solution with fi¬ 
nite (s) > 0. Only the latter, however, is stable, as is well 
known from the mean field theory of the Ising model [24j . 
In the low temperature regime /3 —> oo (with ( k ) (s) > 0) 
the stable solution becomes (s) = 1. In this limit, the 
Fermi function approaches a Heaviside function at x = 0, 
always crossing the linear function g(x) when x —> 0 _ , 
which amounts to ( k) —>• (fc) m i n = 2 d. This is the topol¬ 
ogy phase of the model, in which the ground state of 
the model is a low-clustering universe graph with effec¬ 
tive Hausdorff dimension d e s = (k)/ 2 > 4. Between the 
two phases there is a ferromagnetic phase transition in 
which space-time and the discrete universe emerge spon¬ 
taneously. The exact position of the phase transition can 
be determined by solving numerically the system of equa¬ 
tions m and ®. For example, for d = 4 it is located 
at T cr = (k) CI = 0.4681V. 

In the topology phase of the model the average con¬ 
nectivity ( k ) of the universe increases with decreasing 
temperature. This corresponds to a topological expan¬ 
sion in which the universe actually “unfolds” rather than 
expands. There is no big bang, space-time emerges collec¬ 
tively in a phase transition as a very tight ball of “hyper- 
connected” points with scaling factor D mm oc ./V 2 /^ cr . 
Because of the very high connectivity at the phase tran¬ 
sition, (k) cr = 0(N) this is a finite quantity in the limit 
N —> oo. Note also that there is no singularity at the 
phase transition, the universe graph emerges with an av¬ 
erage distance of 0(1). 

For very high and very low temperatures it is possible 
to derive how the average distance on the universe graph 
scales with the temperature. For very high temperatures, 
P —> 0 one can expand the exponential in the Fermi func¬ 
tion, for very low temperatures /3 —> oo, instead, one can 
replace the Fermi function (1321) with exp(a;) (bear in mind 
that a; is a negative quantity). This gives 

(D) OC Drain N^r , T ->• oo , (33) 

(D) oc D max N~ aT , T -> 0 , (34) 

with a = (4 J 2 /(J + l) 2 ) log (2 J(N — 1)/J + 1) and 

-Dmax ex N 1 ^ for d > 4 and D max oc IV 1 / 4 for d < 4. The 
important point is that, during topological expansion, 
the average distances in the universe scale exponentially 
with the temperature (or its inverse). A tiny decrease 
in temperature causes a large increase in the “topologi¬ 
cal scale factor” of the universe graph. This starts out 
as an almost completely connected graph and expands 
exponentially (with temperature) to a regular graph of 
dimension 4 (if d < 4). 

Let me finally mention a very important point. Mean 
field theory addresses only the average degree < k > of 
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the universe graph. In the phase of interest, however, 
< k ><^C N and thus the exact solution will assign non¬ 
vanishing probabilities to only approximately < k > of 
the possible N(N — l)/2 edges at every vertex. What is 
more, when < k > 7 ^ integer these non-vanishing proba¬ 
bilities will typically differ from each other, a situation 
most pronounced for < k >= 2d + e, where I expect 
some edges with probabilities much smaller than others. 
Finite temperature is thus a mechanism leading to pos¬ 
sible anisotropies in the universe graph. 


V. TOPOLOGICAL BLACK HOLES 

Until now I have concentrated on emerging universe 
graphs with homogeneous space-time “magnetization” 
and vertex degrees, both at T = 0 and at finite tem¬ 
perature in the mean field limit ( k ) ^>1. In this section, 
instead, I will explore possibly inhomogeneous configura¬ 
tions that can emerge at generic, lower values of ( k ). 

To do so I will consider a Peierls droplet of n -C N 
wrong space-time Sj = — 1 spins in an otherwise homoge¬ 
nous universe of N vertices. In the usual Ising model, the 
vertex connectivity is fixed and, when it is high enough, 
these droplets disappear below a finite critical tempera¬ 
ture since the entropy of the droplet becomes lower than 
its internal energy, so that the free energy is dominated 
by the latter [37|. This is given by the boundary con¬ 
tribution E = Eh + 6 where Eh is the energy of the ho¬ 
mogeneous universe graph and 5 is the number of links 
between spins of opposite signs. The internal energy is 
thus minimized by turning all wrong spins so that 5 —> 0. 
This is the standard Peierls argument to explain why all 
the spins align when lowering the temperature and it is 
exactly what happens in the homogenous mean field so¬ 
lution discussed in the previous section. 

In the present model, instead, the links are themselves 
dynamical variables and there is another way to lower 
the internal energy of the droplet. Indeed, turning one 
single spin does not lower the internal energy but it actu¬ 
ally increases it by A E = (k) — 2, where (k) is the mean 
connectivity at the considered temperature, since for one 
’’corrected” link there are now, on average (k) — 1 new 
wrong ones. For high enough temperatures, (k) > 2 and 
the energy cost is a positive quantity. Of course, if the 
links are given external parameters, the only way to lower 
the energy is a simultaneous turnaround of all the wrong 
spins. Here, however the system could simply severe 
the links between spins of opposite signs, thereby lower¬ 
ing the internal energy by A E = — S/2 (1 + J((k) — 1)). 
The resulting configuration is a universe represented by a 
graph containing a hole within which there is simply no 
space-time. This hole is represented by a disconnected 
subgraph for which all vertex spins have value s* = — 1. 
The original universe graph contains a “Schwarzschild- 
like” boundary on which the spectral dimension is lower 
by one than in the rest of the universe graph: this bound¬ 
ary is constituted by all Sj = +1 vertices for which a link 
to a Si = — 1 vertex inside the hole has been severed. I 


will call such a configuration a topological black hole. It 
is characterized by the absence of space-time inside the 
hole, by a lowered spectral dimension on its boundary 
and, at least for integer d , by an entropy proportional to 
the “area”, i.e. the number of vertices on the boundary. 
For d = 2 this is the original Peierls argument [24| . The 
Peierls argument, however, can be generalized to any in¬ 
teger d > 2 37 ], the entropy always scales as the number 
of vertices on the higher-dimensional contour around the 
hole, i.e. its “area”. 

The question arises if a topological black hole is, on 
average, a stable configuration for temperatures down to 
T = 0. The only way to detect the presence of the topo¬ 
logical black hole from within the main universe graph is 
through the lower spectral dimension on its boundary. 
And indeed, the lowest-energy instabilities come from 
vertices on the boundary forming links between them¬ 
selves. Forming such a link would lower the system en¬ 
ergy by 1/2 via the ferromagnetic vertex term in the en¬ 
ergy, on the other side, it would cost energy due to the 
antiferromagnetic link term. In total, one such link would 
cost an energy 


A E = J((k) - 1) - 


1 


(35) 


Note that forming links to other vertices not on the 
boundary costs more energy since they possess already 
more links. Topological black holes are thus, on average, 
stable configurations against radius connections if (l35l) is 
a positive quantity. Using J = l/(4d — 1) this stability 
condition can be formulated as 

{k)>2d+\=k{T = 0) + \. (36) 


This shows that topological black holes are typically sta¬ 
ble down to very low temperatures. 


VI. FUTURE DIRECTIONS 

Two further steps are needed to make definitive con¬ 
tact with gravity. First and foremost the renormalization 
flow of the dimension and the local limit at the IR fixed 
point have to be derived and it must be shown how the 
Einstein-Hilbert action arises there. Secondly it remains 
to show how the local limit admits a causal structure and 
why the emerging continuous space is Lorentzian. These 
are the subjects of further investigation. 
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